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Abstract 

The electroweak baryogenesis depends on the profile of the CP-violating bubble 
wall created at the first order phase transition. We attempt to determine it by 
solving the coupled equations of motion for the moduli and phases of the two Higgs 
doublets at the transition temperature. A variety of CP-violating bubble walls are 
classified by boundary conditions. We point out that a sufficiently small explicit 
CP violation gives nonperturbative effects to yield the baryon asymmetry of the 
universe. 
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1 Introduction 



To generate the baryon asymmetry of the universe (BAU) starting from a baryon- 
symmetric universe, Sakharov's three conditions must be met. If there was no B — L 
down to the electroweak era, baryogenesis mechanism should work at the electroweak 
phase transition (EWPT) for the present universe to be reahzed|jl[]. Even if some B — L 
existed, what happened at the EWPT would affect the present BAU. Among the condi- 
tions, to have sufficient CP violation necessary for the BAU will require some extension 
of the Higgs sector in the minimal standard model. Further, the minimal standard model 
satisfying the present lower bound on the Higgs mass cannot make the EWPT first order 
and cannot suppress the sphaleron process after that. Two-Higgs-doublet extension of the 
standard model is the simplest extension, which can yield the Higgs-sector CP violation, 
and includes the Higgs sector of the minimal supersymmetric standard model (MSSM) as 
a special case. Here we concentrate on such two-Higgs-doublet models. 

At the first-order EWPT, a bubble wall dividing the broken and symmetric phases is 
created and grows to convert the whole universe into the broken phase. The profile of the 
bubble wall will be characterized by the classical configuration of the Higgs and gauge 
fields. It is essential in any scenario of the electroweak baryogenesis to know the spacetime- 
varying profile of the bubble wall, in which the relative phase of the Higgs scalars violates 
CP. In literatures, some functional forms of the profile were assumed00]. They should 
be determined, however, by the dynamics of the gauge- Higgs system. That is, the bubble 
wall will be a solution to the classical equations of motion which has the effective potential 
as the potential for the Higgs sector. The observable CP violations in our world offer the 
boundary condition in the broken phase for these equations of motion. 

In a previous paper |Q, assuming some appropriate form of the effective potential, we 
found a solution such that CP-violating relative phase becomes as large as 0(1) around 
the wall while it completely vanishes in the broken and symmetric phase limits. In the 
solution, the moduli were fixed to be the kink shape. Such a solution, which we call 
as a "solution with kink ansatz", is interesting since it does yield sufficient amount of 
the chiral charge flux through the wall surface, which will be turned into the baryon 
number in the symmetric phase region by the sphaleron transition. The solution emerges 
when a mechanism similar to the spontaneous CP violation operates in the intermediate 
region. For such a mechanism to work, some of the parameters in the effective potential 
are restricted similarly to the case of the spontaneous CP violation at zero temperature, 
although less constrained at finite temperature. 

The next task will be to determine the moduli themselves of the Higgs fields dynami- 
cally as well as the relative phase of them. In this paper, we solve the coupled equations 
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of motion for the moduli and the relative phase of the two Higgs fields by imposing a 
discrete symmetry between the two Higgs doublets to reduce the dynamical degrees of 
freedom. With the same parameters of the effective potential as those with the kink 
ansatz, we present a solution, which we call as a "solution without kink ansatz" , having 
larger CP- violating phase and non-kink type moduli. As expected, this solution has lower 
energy than that with the kink ansatz. 

We also give another class of new solutions, which are caused by a similar mechanism 
of CP violation in the intermediate region as above but are free from the constraints of 
zero-temperature spontaneous CP violation. The nontrivial phase dependence of cubic 
terms in the moduli, which is essential to the first-order EWPT, plays a crucial role in 
this mechanism. Such terms will be induced from bosonic loops, especially those of the 
scalar partners of quarks and leptons in the MSSM and the extra Higgs scalars. 

If CP violation is completely spontaneous, a CP-violating bubble is created with the 
same probability as its CP-conjugate one. Then the net baryon number would result in 
zero. We show that a tiny explicit CP violation, which is consistent with the present 
bound on the neutron electric dipole moment (EDM), can nonperturbatively resolve the 
degeneracy between the CP-conjugate pair of the bubbles to leave a sufficient BAU after 
the EWPT. 

In section 2, we postulate the effective potential of the Higgs sector and present the 
equations of motion. In section 3, we classify possible solutions according to the boundary 
conditions, and give some examples of "solutions without kink ansatz". We also investi- 
gate solutions in the presence of explicit CP breaking, and point out the importance of 
the breaking to generate the BAU. Section 4 is devoted to conclusion and discussions. 

2 The Model 

One of the most attractive features of electroweak baryogenesis is the fact that it depends 
only on physics which can be tested by experiments at present as well as in the near 
future. We emphasize that it is CP violation at the first-order EWPT, not that observed 
in laboratories, which affects the generated baryon asymmetry. The latter is restricted 
by various observables such as the neutron EDM, while the former is not directly. Here 
we attempt to relate the observable CP violation to that at the EWPT. 

2.1 Ansatz for the effective potential 

Given a model, one can calculate the effective potential near the phase transition tem- 
perature, which determines the order of the EWPT as well as its dynamics. We expect 



3 



that the bubble wall will be well approximated by a solution to the classical equations of 
motion derived from the effective Higgs potential, which consists of the tree-level poten- 
tial with radiative corrections and finite-temperature contributions. Even if the tree-level 
Higgs potential is CP conserving, the corrected potential parameters may predict sponta- 
neous CP violation, or the corrections may bring effects of various explicit CP- violating 
parameters, such as those in the soft supersymmetry breaking parameters, into the ef- 
fective potential. The CP violation in the Kobayashi-Maskawa matrix could also induce 
CP- violating term in the effective potential through higher order diagrams. 

Here we postulate a functional form of the effective potential, which is equipped with 
the feature of first-order EWPT. We assume that it is well approximated by a gauge- 
invariant polynomial up to the fourth order of the moduli of the Higgs fields. The most 
general gauge-invariant Higgs potential is 

+ ^Ai($l<l>i)2 + ^A2($5$2)' + A3($I*l)(*2$2) (2.1) 
- A4($I$2)($5$l) - |iA5($I$2)' + [A6($I$l) + A7(<I>5$2)]($I$2) + h.C.j , 

where , m^, Ai, A2, A3, A4 G R and m|, A5, Ae, A7 G C, three of their phases are in- 
dependent and yield the explicit CP violation. To avoid the tree- level Higgs- mediated 
flavor-changing neutral current (FCNC) interactions naturally, one may impose the dis- 
crete symmetry on the lagrangian[0, which is broken by Aej 7^ and softly by ml 7^ 0. 
The Yukawa interactions should have that symmetry not to induce divergent Aej-terms. 
In the MSSM, the parameters in the Higgs potential are given by 

Ai = A2 = \{9' + g'\ A3 = - A4 = \g\ 

A5 = Ae = A7 = 0, (2.2) 

where g {g') is the SU{2) {U{1)) gauge coupling, the mass parameters are completely 
arbitrary, and the Yukawa interactions respect the discrete symmetry. Hence CP is not 
violated by the Higgs sector at the tree level in the MSSM. Nonzero m\ will induce finite 
corrections to As^ej, which violates CP explicitly or spontaneously. 

At finite temperature, there arise new kind of corrections, which lead to the first-order 
EWPT. In general, the one-loop finite-temperature corrections to the effective potential 
have a form ^ 

V{v-T) = ^^^dAlA{a\), (2.3) 

where A stands for particle species running through the loop, 

Ib,f{(^^) = dxx^ log (2.4) 
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with ttA = fnA{v)/T, v being the VEV of the Higgs fields, and counts the degrees of 
freedom with sign (e.g., dw = Q, dz = 3 and dt = —6). One can expand Ib,f{,o?') around 

= to observe that the expansion of Ib contains (a^)'^/^-term, which originates from 
the Matsubara zero mode of the loop momentum^]. Hence if m?{v) ~ at f ~ 0, the 
effective potential has terms with the correct sign, which is the characteristic of the 
first-order phase transition. Here we assume that the effective potential at the EWPT 
is approximated by the polynomial of the form of (|2.1|) plus t>^-terms irrespective of the 
validity of the perturbation and the high-temperature expansion. 

First of all, we shall concentrate on the case with no explicit CP violation in the 
effective potential. Then the ansatz for it is 

Kff(pi, 0) = ]^m\pl + ^mlpl + mlpip2 cos 6 



Ai 



X 



A3 — A4 — A5 cos(26') 2 2 
A P1P2 



-lj^i.^fiPlP2 + A7P1P2) COS 6* 

'Ap\ + (5o + Bi cos 9 + B2 cos(20))p2p2 

+{Co + Ci cose + C2Cos{29))pipl + Dpi 



(2.5) 



where we parameterize the Higgs VEV asQ 



6 — 61 — 62- 



(2.6) 



As we noted, the parameters in (|2.5| ) should be understood to contain radiative and finite- 
temperature corrections, which can be calculable once the model is specified. According 
to the high-temperature expansion, Ba and Ca may arise from the bosonic loops, among 
which ^-dependent terms, B12 and Ci 2, could arise only from those of squarks, sleptons 
and Higgs scalars, that mix scalars and pseudoscalars in the MSSM. Viewed as a function 
of cos^, the effective potential is written as 



V,s{p.,9) 



A 



''plpl + 2{B2pip2 + C2Plpi) 



X < COS 61 + 



+ 



-2ml + Aepf + Arpj + 2{Bipi + C1P2) 
2A5P1P2 + 8(52Pl + C2P2) 



(2.7) 



When 



A5PIP2 + 4(S2ptp2 + C2P1P2) < 0, 



^Here we assume that /7(l)em gauge symmetry is unbroken at the EWPT. 



(2. 
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and 



-2ml + ^6P? + ^ipI + 2(5ipi + C1P2) 



< 1, 



2A5P1P2 + 8(52pi + C2P2) 
CP is spontaneously broken. At zero temperature, i?a = = and (pi,P2) 
f sin /?o) so that these conditions are reduced to the well-known ones, 

-2m\/vl + Ae cos^ /So + A7 sin^ /?o 



A5 < 0, and 



2A5 sin I3q cos I3q 



< 1. 



(2.9) 

{vo cos/3o, 

(2.10) 



Here all the parameters are evaluated at zero temperature. The former condition (A5 < 0) 
is hardly fulfilled, as long as A5 > at the tree level, since almost all radiative correc- 
tions except for those by the gauginos and Higgsinos contribute to make A5 positive in 
the MSSM[0]. Even if these conditions are satisfied, such a model will be bothered by 
a light scalar, which is inevitable in the presence of spontaneous violation of a discrete 
symmetry^]. In contrast, the conditions ( |2.8|) and ( p.9|) can be satisfied in the intermedi- 
ate range, since pi (P2) varies from v cos (3 [v sin (3) to zero as going from the broken phase 
to the symmetric phase region. Since the finite-temperature corrections both of bosons 
and fermions to A5 are positive at the one-loop level, requiring ( |2.8|) and ( |2.9| ) to be satis- 
fied only near the bubble wall and not to be satisfied in the broken phase would constrain 
some parameters in the model. Especially in the MSSM, for the finite-temperature cor- 
rections to negative A5 to be small, some soft-supersymmetry-breaking mass parameters 
will be bounded from belowQ. 

If CP is violated spontaneously, there are two kinds of bubbles {pi,0) and {pi,—6) 
which are CP conjugate with each other and have the same energy density. Such bubbles 
would be nucleated with the same probability, so that no net baryon number would be 
generated. An explicit CP violation will resolve this degeneracy and leave the finite BAU. 
Although various types of explicit CP violation in the effective potential are possible, we 
restrict to a simple case that only the mg is complex and the other parameters are real. 



Hence we replace m\pip2 cos 6 in Vcs of ( p.5|) by 

mlpip2 cos{6 + 6), 



{ml e R) 



(2.11) 



In the MSSM, nonzero 6 is induced from the CP violation in the scalar trilinear terms 
and \6\ < 10^^ does not contradict with the present bound on the neutron EDM||TO| . 



The effect of 6 on the baryon asymmetry was discussed in [|T^] perturbatively. We shall 
see later that 6 nonperturbatively discriminates the CP-conjugate pair of the bubbles for 
some types of solutions. 
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2.2 Equations of motion 

We suppose that the dynamics of the bubble wall is governed by the classical equations 
of motion derived from the lagrangian 

£ = -^F^^F^f^^ - -B^^B^'' + {D,<^d'' D^<^, - $2; Tc), (2.12) 

^ ^ i=l,2 



where, 



If the bubble is spherically symmetric and sufficiently macroscopic, it is approximated 
as a planar object so that the system is reduced to one-dimensional. Further, when the 
bubble grows with a constant velocity keeping the shape of the critical bubble, it can 
be viewed as a static object. For this effective one- dimensional system, the gauge fields 
may be written in pure gauge forms and a solution in this case will have lower energy 
than solutions with nontrivial gauge configuration. Hence we assume that the bubble 
wall consists of only the Higgs scalar, and work in the gauge where the gauge fields are 
eliminated. Since the fermions couple to gauge-noninvariant phases of the Higgs instead 
of the invariant 6*, the CP violation in the Yukawa interactions is made unambiguous once 
the gauge is fixed in this manner. 

Then the bubble wall profile is characterized by the space-varying order parameters 
{pi{z),6i{z))^ which satisfy the following equations of motion, 

where z is the coordinate perpendicular to the wall. From the requirement that the gauge 
fields are pure-gauge type, a constraint equation, "sourcelessness condition", should be 
fulfilled: 

P?(=)^^+P^W^^ = 0, (2.15) 

dz dz 

Now our task is to solve the equations ( p.l3| ) and ( p.l4| ) with the constraint ( |2.15| ) satis- 
fying certain boundary conditions in the broken phase and symmetric phase regions. 

Before solving the equations, we reduce the dynamical degrees of freedom by adopting 
some ansatz. First, we require that the moduli of the Higgs scalars take kink shapes of 
the common width when 9 is sufficiently small. In practice, we assume that when 5 = 0, 
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9{z) 



V COS (3 

V sin (3- 



1 + tanh(a2;) 
2 

tanh(az) 



(2.16) 



are solutions to (|2.13|) , where 1/a is the common wall thickness. This restricts some of 
the parameters in the effective potential ( p.5|) ||^. Then the equation of motion for d{z) in 
the background of the kink Pi{z) is reduced to 



y {^-y) 



+ y{l-y){l-Ay) 



de{y) 



dy"^ ' dy 
= 6sin(0(y) + 5) + [c(l -y f-e{l- y)] sme{y) 

+ [^(l-y)2-2/(l-y)]sin(2%)), 

where we introduced a dimensionless finite-range parameter y by 

1 



y 



[1 — tanh(a2;)) 



and the parameters in ( p.l7| ) are defined by 



Aa? sin (3 cos (3 ' 



V 



c = 



32a2 

,2 



(Ai cot^ (3 + \2 tan^ (3 + 2(A3 - A5)) - 



2 sin^ (3 cos2 (3 



V 



d = 



V 



■(Ae cot /3 + A7tan/3), 



/ 

e + f 



Aa^ ysin/? cos/5^ 
' + 



Isin/? cos/? 



; {A cos^ /5 + 5o cos^ (3 sin /3 



4a2 



sin^ f3 cos^ /? ^ 
+ Cocos/?sin^/? + Dsin^/? - — ) 

V 

with A3 = A3 — A4. Here 6{z) is related to ^1,2(2:) by 

ei{z) / sin^ (3 = -e2{z) / cos^ (3 



(2.17) 



(2.18) 



:2.19) 



(2.20) 



which means that 0i{z) and 02{z) satisfy the condition ( |2.15| ) when Pi{z) is the kink 
shape ( p. 16 ). We refer the solution to ( p. 17 ) as "the solution with kink ansatz". We 
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found several solutions for various boundary conditions when / = in Among them, 
when d < (A5 < 0), there is an interesting solution which connects the CP-conserving 
vacua and largely violates CP in the intermediate range around pi{y) ~ 1/2 by satisfying 
( p.8|) and (|2.9| ). The solution has so large a phase as 9{y = 1/2) > 0.3 to generate 
sufficient chiral charge flux. For the same set of parameters, we also found that a tiny 



explicit CP violation of the form ( p.ll|) can nonperturbatively enhance the lower-energy 



bubble over its CP-conjugate partner to leave sufficient BAU [|T2| 



The large 6 means, however, that the kink ansatz may no longer be adequate. Thus 
we have to solve coupled equations for the moduli and the phases. In order to simplify 
the analysis, we impose a discrete symmetry on the effective potential under 

Pi ^ P2, ^1 ^ -02, (2.21) 

and assume that it is not spontaneously broken to reduce the dynamical degrees of free- 
dom. This means that the parameters in (|2.5|) should satisfy 



ml = ml = m^, Ai = A2, Ag = Ay, 

A = D, Bq = Co, Bi = Ci, B2 = C2, 
and that tan /3 = 1 by definition. For 

Piiz) = P2{z) = p{z)/V2, ei{z) = -e2{z) = e{z)/2, (2.22) 

we have two coupled equations: 



d'p{z) 1 , , (de{z)V dV, 



while the "sourcelessness condition" is automatically satisfied. Here the effective potential 
is written in terms of the reduced numbers of the parameters as 



1 1 ~ 

VeE{p:0) = - (m"^ + mlcos{e + 6)) p'^ -\ [A - A5 cos(26') - 4A6 cos6' 

2 ^ ^ 16 L 



— ^[A + Bo + Bicos9 + B2 cos(2e)] p^ (2.25) 
v2 

where A = Ai + A3 — A4. The parameters in ( |2.25| ) will be further constrained by requiring 
that Vcs should have two degenerate minima corresponding to the broken and symmet- 
ric phases. We relate the parameters to those in the potential with the kink ansatz. 
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{b,c,d,e, f), to compare "the solutions without kink ansatz" to those with it. That is, 
some of the parameters in ( p.25| ) are related to (6, c, d, e, /) directly by 



A 



A5 + ^^(c + 2), 



4a2 , 

V 



Bo 



/, 



V V 

and from the kink ansatz for the case of ^ = 0, the others are given by 



4a' 



m 



3; 



V 



Bi — B2, 

2 4a2 



-c. 



(2.26) 



,v 2 2- (2-27) 

Once {b,c,d,e, f), the wall thickness 1/a and the Higgs VEV v at the transition tem- 
perature are given, the effective potential is determined. Such an effective potential has 
the degenerate minima and covers a restricted class of the potentials of the form ( |2.25| ). 
We expect that the solutions obtained below will survive when we slightly deform the 
effective potential from the restricted class. 



2.3 Energy density and BAU 

The nucleation rate of the bubble is roughly given by expl—AnR^S/Tc), where Rq is the 



radius of the critical bubble, which is approximated by Rc — y 3Fc/(47raf 2). Here Fc is 
the free energy of the critical bubble and is found to be about (140 ~ 160)T|l^. S is the 
energy density of the bubble, 

2 / \ 2" 



For "the solution with kink ansatz" , this is given by 
,2 



+ Ves{puP2, 0) 



^2.28) 



^ " 3 



+at>2 sin^ (3 cos' (3 I dy 
Jo 



dO 



y{i-yY\-r\ + ^6 [1 - cos(^ + 5)] 



dy 



y 

+ [^(l-y)-/][l-cos(2^)]^ 



(2.29) 
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where is the energy density of the trivial kink solution. It was evaluated for several 

solutions in and With the discrete symmetry, the energy density is reduced to 




£ = dyl ay{l - y) 



In the charge transport scenario]^, the generated baryon number is proportional 
to the hypercharge flux Fy, which is caused by the difference of the reflection rates of 
chiral fermions by the CP- violating bubble wall. That is, 

'^^N^-nat,-^,-rT, (2.31) 

where s is the entropy density, M a model-dependent constant of 0(1) and k the coefficient 
of the sphaleron transition rate in the symmetric phase given by k = 1.09 ± 0.04[Q. r 
may be approximated by the mean free time or diffusion time D/vw with the diffusion 
constant D of the charge carrier, where D ~ 1/T for quarks and D ~ (10^ ~ 10'^)/T 
for leptons[|T5|. It was shown that if the forward scattering is enhanced, even for the 
top quark, tT ~ 10 ~ 10^ depending on f^,, where the maximum value is realized at 
Vu, — l/-\/3 10. For this optimal case, 

^ - 10-' ■ (2-32) 

Then the hypercharge flux Fy/{vJT^) > 0(10^^) would be necessary to explain the 
present BAU. 

As we noted, if CP is violated spontaneously, the bubble with the profile (pi, 9) degen- 
erates with its CP-conjugate partner (p^, —9). Then their nucleation probability and dy- 
namics will be completely the same so that the net baryon number generated is completely 
canceled with each other. Even if an explicit CP violation will resolve the degeneracy, 
two (or more) kinds of bubbles contribute to the baryon number in opposite signs. When 
we denote the baryon number generated by the j-th kind of bubble as {nB/s)j and the 
nucleation rate of it as Nj = exp(— 47rP^£^j/Tc), the BAU may be approximately given 
by 

where the expansion rate of each bubble is assumed to be almost the same. We shall 
estimate it for some numerical solutions obtained in the next section. 
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3 Classification of Bubble Walls and Numerical 
Examples 



We show some solutions to the coupled equations ( p.23|) and ( p.24| ), that is, "the solutions 



without kink ansatz". First, one must specify the boundary conditions in the broken and 
symmetric phases. 

By construction, V"cfr(p, 0) has two degenerate minima. In the broken phase, the bound- 
ary value 6*0 in (f , 6*0) = (p(y = 0),6{y = 0)) is determined by minimizing V^fj(t>, 6*0). 
Namely, for a fixed v, 6q is the solution to 

^(As cos 6*0 + \e)v^ + \/2v{Bi + AB2 cos 6*0) 



sin (6*0 + 5) 



sin^o = 0. (3.1) 



2 

When the parameters are written in terms of {b,c,d,e, f) by (|2.26|) and ( p.27| ), this is 
expressed as 

b sin(6'o + 5) + (c + cos 6*0 - e - 4/ cos 6*0) sin 6*0 = 0. (3.2) 

In general, for small \6\, this has a solution 6q = 0{6) (mod tt). This corresponds to the 
CP-conserving minimum when 5 = 0. For 6 = and d—Af < 0, if \{b + c — e)/{d — 4/) | < 
1, there exists a nontrivial minimum cos 6*0 = —{b + c — e)/{d — 4/), which implies spon- 
taneous CP violation in the broken phase. Such a boundary value was considered in our 
previous workQ under the kink ansatz. Even if 6*0 for 6 = obtained above is very small, 
for example, O(10~^), a small 6 happens to drive 6q nonperturbatively to a large value, 
such as 0(10^^), as seen by solving ( [3.21 ). If the parameters in the potential match this 
condition only at finite temperature, this will be a new possibility to have sufficient CP 
violation at the EWPT. 

As for the symmetric phase, it is obvious that p{y = 1) = 0. Although 61 = 6{y = 1) 
seems to be undetermined, it is not the case when mg 7^ (6 7^ 0). This is because the 
equation for 6{y), (|2.24|) written in terms of y, has singularity at y = 1, which makes the 
solution singular there unless sin(6'i + 6) multiplying the most singular m|-term vanishes. 
For such a singular solution the energy density diverges so that it cannot be realized 
dynamically.^ Hence the boundary condition for 6{y) in the symmetric phase is 

9^ + 6 = mi, (nGZ). (3.3) 

As long as the effective potential is parameterized as in the previous section and \6\ 

is sufficiently small {\6\ < 10^'^), the possible lowest energy solutions will be classified in 

terms of the boundary values of 9{y) as in the table below. 

^Strictly speaking, since the kink ansatz is imposed when 0{y) = 0, the effective potential satisfying 
the ansatz has a minimum at a different point from the prescribed one when 6^0. 



3t 



In practice, the integration region in ( 2.28 ) is finite so that the energy density is finite. We, however, 



expect such a solution is hard to be nucleated because of its large energy. 
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Re m -Im m 

■-1.0x10"^ 

-8.0x10"" 

-6.0x10"" 

-4.0x10"" 

-2.0x10"" 



0.0x10" 

-4-3-2-101234 

x=a7. 

Figure 1: The wall profile for (6, c, (i, e, /) = (3, 5, 5, 7, —1.25) and 5 = 10^^ as a function 
of the dimensionless variable x = az. Note that the CP-violating imaginary part ( dotted 
curve ) is much smaller than the real part ( solid curve ). 






6 = 


6^0 


(A) 


9{y) = (trivial solution) 


e{y) = 0{6) for 


(B) 


= ^1 = 

spontaneous CP violation in the wall 


00 = 0{6) and Oi = -6 
\9{y)\ » \5\ at y ~ 1/2 


(C) 


00^0 and = 
spontaneous CP violation in the broken phase 


l^ol » 0{6) and 6^ = -5 
\eiy)\ » \5\ aty< 1/2 


(D) 


= and 6i = tt 
maximal CP violation in the wall 


9o = 0{6) and ^1 = TT - 5 
maximal CP violation in the wall 



As an example of type (A) with 5 7^ 0, we study the case for (6, c, d, e, /) = (3, 5, 5, 7, 
— 1.25). We plot the wall profile or the mass function in the Dirac equation, (psin(6'/2), 
pcos{6/2)), in Fig. 1 and the chiral charge flux obtained from it in the dimensionless unit 
of Fq/{vwT^ {Ql — Qr)) in Fig. 2, where the wall velocity and temperature are taken to 
be f ^ = 1/ -\/3 and T = lOOGeV, respectively. Q denotes a certain quantum number such 
as the hypercharge. For this type of solution with Oiy) = 0(5), it seems difficult to have 
sufficient BAU except for a restricted set of the carrier mass and the wall thickness. 

If a solution of type (D) exists, it will incorporate so large CP violation near the 
bubble wall to generate sufficient BAU. We have attempted to find such a solution by 
adjusting the parameters in the effective potential, but have not succeeded. This may be 
because the cost of the derivative term in the energy cannot be compensated by the gain 
of the potential term along such a path from the broken to the symmetric phase. 

If type (D) solution does not exist, the only possible one accompanying large CP 
violation 9 ~ 0(1) in an intermediate range will be of the type (B) or (C). The range 
where 9 is large should not be near the symmetric phase, where p{y) ~ 0, since the CP 
violation in the Yukawa interactions is proportional to piy) sui(9{y) /2). A solution of type 
(B) for 5^0 was found under the the kink ansatz [ll2i . 



13 



log (nvT) 



2 



-1.5 



-0.5 



0.5 



1.5 



-2 



-1- 







1- 



-3 




7 



logjo(a/T) 



Figure 2: Contour plot of the net chiral charge flux, normalized as 
logj^Q [—Fq/{vwT'^{Ql — Qr))] for the profile shown in Fig. 1. Here we take Vyj = 0.58 
and T = lOOGeV. 

Let us denote the solutions of the CP-conjugate pair in the case of 5 = as (p, 6) and 
{p,—9). For sufficiently small \6\ ^ 0, there generally appear three solutions, (p+,^^+), 
{p^,9^) and (p°,6'°), such that they respectively approach {p,0), {p,—9) and the trivial 
kink solution (pkink, 6* = 0) as 5 ^ 0. Note that (p", 9^) has no CP-conjugate partner. 
An example presented here employs the same parameter set and the boundary condition 
as in |12| to solve the equations ( |2.23| ) and ( p.24| ) for various 6 starting from 6 = 0. 
That is, (6, c, d, e, /) = (3, 12.2, -2, 12.2, 0) and 9i = -6 while 9o is determined by Q. 
There are three numerical solutions, {p^,9'^), {p^,9^) and {p^,9^), in the order of the 
decreasing energy, satisfying the same boundary conditions for 6 < 0.0045. For larger 
6, only {p^,9^) exists. Their energy densities, measured from that of the trivial kink 
configuration, 6S = S — av'^/3, are plotted in Fig. 3 in the unit of av"^. To illustrate 
how CP is violated in the intermediate region, we show the contour plot of the effective 
potential together with the path {p~^,9~^) for 6 = 10~^ in Fig. 4. The mass functions of 
each profile for 6 = 10~^ are plotted in Fig. 5. For illustration to evaluate BAU, we take 
the case of 5 = 10^'^. Since the energy densities of them are. 



SS[p^9'] 
SS[p-',9^] 
6£[p-,9-] 



1.21567 X lO^W, 
-5.6350 X 10~^av^ 



(3.4) 



1.2012 X lO^W 
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Figure 3: Energy densities of the bubble walls for (6, c, d, e, f) = (3, 12.2, —2, 12.2, 0) 
measured from the trivial kink solution in the unit of av^. 




Figure 4: Contour plot of the effective potential together with the path of {p^,9^) con- 
necting the two vacua for the same parameter set as Fig. 3 and 6 = 10^^. 




x=az 
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Figure 5: The wall profiles or the mass functions for (6, c, d, e, /) 
and 6 = 10"'^ as functions of x = az. 
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Figure 6: Contour plot of the net cliiral charge flux, normalized as 
logio [~Fq/ {vyjT^iQL — Qr))] for the profile shown in Fig. 5. Here we take = 0.58 
and T = lOOGeV. 



the ratios of the nucleation rates are 

-± = 12.23, — = 196.0, (3.5) 

iVo iVo 

where we have used the expression for the nucleation rate of the critical bubble 

A'.exp(-l^). (3.6) 



with Rc — y3Fc/{4:Tcav'^) and Fc — 145Tc. By neglecting an extremely small contri- 
bution from {p^,6^), as understood from the small imaginary part in Fig. 5, it is found 
that about 88 % of baryon number generated by (p~,0~) survives the cancellation by 
(p"^, 6^). Thus 6 = 10"^ will be sufficient to leave the BAU for wide range of the carrier 
mass and the wall thickness. As seen from Fig. 3, only one species of the bubbles exist 
for 6 > 0.0045 so that the cancellation among the bubbles does not occur. The net chiral 
charge flux ( |2.33| ) for 6 = 10^^ is presented in Fig. 6. 

Some solutions of type (C) for 6 = were found under the the kink ansatz[Q. The 
profile for non- vanishing 6 is similar to those presented there [Q, which monotonously 
connects {v,6q) to (0, = —5). The magnitude of the chiral charge flux is controlled by 
that of 9n. 
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4 Conclusion and Discussions 



In conclusion, assuming that the EWPT is of first order accompanying nucleation and 
constant growth of the broken-phase bubbles within the symmetric phase and that the 
dynamics of the bubbles are governed by the classical equations of motion of the gauge- 
Higgs sector, we have given some examples of the bubble profile by solving the equations 
of motion derived from the effective potential. According to the boundary conditions, the 
possible types of solutions are classified. Among them, the profiles of type (B), (C) and 
(D) will yield large CP violation around the bubble wall, though the type (D) seems hard 
to be realized dynamically. The solutions of type (B) require the parameters to admit 
spontaneous CP violation in some region between the broken and symmetric phases. 
The explicit CP violation, which is small enough not in contradiction with observations, 
can avoid the cancellation of the generated baryon number between the CP-conjugate 
bubbles, leading to the BAU. The type (A) bubbles will yield insufficient BAU except for 
a restricted set of the carrier mass and the wall thickness. 

If the temperature-dependent parameters are calculated in the MSSM at the one- 
loop level, the parameter / will be generated only by higher order terms in the high- 
temperature expansion of the squark or slepton contribution and the condition dp^ — 
4/p < will be satisfied by c/ < 0. For d = f^/(4a^) ■ A5 < to be realized, we must have 
tree-level negative A5 in the two-Higgs-doublet model or radiatively induced negative A5 
in the MSSM. In this case, we must avoid the problems peculiar to the spontaneous CP 
violation, that is, too light scalar and the domain wall. Both will be cured by a small 
explicit CP violation. In the MSSM, it will be unnatural to assume that all the soft 
supersymmetry breaking parameters do not introduce any explicit CP violation at the 
electroweak scale. As we noted in the previous section, if CP is spontaneously violated 
in the broken phase when 5 = 0, a small 6 nonperturbatively happens to drive large 
CP violation. This will further restrict some set of the parameters in the model. These 
constraints will be examined in a future publication. 

We have presented the chiral charge flux for the profile of type (A) and (B), which is 
calculated supposing that one of the Higgs doublets couple to a fermion mass eigenstate. 
In practice, the weak hypercharge plays its role and the hypercharge flux is additive for 
each fermion species in the standard model. This is because the up-type quark with 
positive hypercharge couples to e*^^ = e*^/^, while the down-type quarks or the leptons 
with negative hypercharge couples to e*^^ = e~*^/^ or (e*^^)* = e~*^/^ in the MSSM or the 
two-Higgs-doublet models with the discrete symmetry to avoid the tree-level FCNC In 
a supersymmetric standard model, there also exist contributions from the superpartners 
with nonzero hypercharge such as the scalar partners of the fermions and those of the 
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Higgs scalars. The effects of the chargino were considered in |]TB 



Finally, we comment on the discrete symmetry (p.21|) , which is a working hypothesis 
to reduce the dynamical degrees of freedom. The numerical method we employed to solve 
the boundary value problem is the relaxation algorithm, which varies the configuration 
to approach the correct solution starting from a prescribed initial configuration. Toward 
our final goal to find solutions to the equations for the full dynamical degrees without the 
discrete symmetry, the knowledge of the solutions obtained here will help us to prepare an 
appropriate initial condition for these equations. That is, as long as the solution with the 
discrete symmetry is stable in the enlarged functional space, we can obtain the desired 
solution by graduately modifying the parameters from those with the discrete symmetry 
to find the next initial configuration, step by step. In fact, we checked that for the solution 
of type (B), the discrete symmetry is not broken spontaneously, that is, the solution with 
the discrete symmetry is still a solution to the equations of motion for all possible order 
parameters. 
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